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ABSTRACT 


Three  integral  equations  are  derived  and  formulated 
for  numerical  solution  for  the  currents  induced  on  a 
resistively  loaded  conical  antenna  over  a ground  plane. 

The  first  integral  equation  is  a relatively  simple  one  for 
a cone  without  a topcap.  The  second  and  third  integral 
equations  are  applicable  to  a cone  with  or  without  a top- 
cap,  but  the  latter  equation  is  relatively  cumbersome, 
involving  complicated  kernels  with  various  singularities. 

A computer  code  has  been  developed  for  each  of  the  three 
methods.  Numerical  data  for  current  distributions,  input 
impedances  and  radiation  patterns  are  presented  for  resis- 
tively loaded  and  unloaded  structures,  both  with  and 


without  topcap 
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SECTION  I 


INTRODUCTION 

In  a companion  report  [ij,  the  static  analysis  of  a 
conical  antenna  over  a ground  plane  is  presented.  In  this 
report,  the  analysis  is  extended  to  treat  the  time-har- 
monic case  and  to  incorporate  a model  of  the  resistive 
loading  of  the  structure.  The  resistive  loading  is  intended 
to  reduce  the  effect  of  diffraction  from  the  cone  edge  at 
the  higher  frequencies. 

If  the  cone  has  no  topcap,  the  analysis  may  be  consi- 
derably simplified  and  a simple  integral  equation  for  this 
situation  is  derived  in  Section  II  and  implementation  of 
a moment  method  solution  is  considered  in  Section  III.  In 
Section  IV,  an  integral  equation  for  the  cone  with  a top- 
cap  is  derived.  Presented  in  Section  V are  numerical  results 
in  the  frequency  domain  currents  for  a loaded  conical  antenna 
both  with  and  without  a topcap.  In  Appendix  A expressions 
are  derived  for  the  computation  of  fields  from  the  currents 
and  Appendix  B gives  the  derivation  of  an  alternate  integral 
equation  from  that  derived  in  Section  IV. 


SECTION  II 


FORMULATION  OF  AN  INTEGRAL  EQUATION  FOR 
A BICONE  WITHOUT  ENDCAPS 


For  the  symmetrically  driven  biconical  structure  of 
Figure  1,  the  current  on  both  the  cone  and  its  image  are 
radially  directed  and  have  no  circumferential  ($)  variation. 
Hence  the  magnetic  field  tangent  to  the  cone  iB  ^-directed 
and  the  boundary  conditions  can  be  satisfied  by  fields  which 
are  transverse  magnetic  (TM)  to  r.  Thus,  the  fields  may 
be  completely  determined  by  a radially-directed  vector 
potential  A - A^?  [2],  In  an  eigenfunction  solution  to 
such  problems,  the  fields  are  determined  in  the  bicone 
region  from  a vector  potential  A^  which  comes  from  a 
homogeneous  solution  of  the  wave  equation.  In  order  to 
derive  an  integral  equation,  however,  A^  must  be  expressed 
in  terms  of  the  current  on  the  bicone.  In  particular,  a 
free  space  Green's  function  is  to  be  found  for  the  vector 
potential  Af  due  to  a unit  r ad ia 1 ly- d ir ec t ed  current  element 
A Buperpoa i t ion  in  egral  then  expresses  the  total  vector 
potential  due  to  currents  on  the  cone. 

Beginning  with  the  assumption  that  the  magnetic  field 
is  determined  from  A - A 
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Figure  1.  Geometry  of  cone  over  a ground  plane. 
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and  using  Maxwell's  equations, 


VxE  -joiyH 


f}- 

f ¥ 


VxH  - jweE+J 


one  readily  determines  by  standard  procedures  that  X satis- 
fies the  vector  Helmholtz  equation 

VxVxA  - k2X  - yJr£  - ju)UeV4>  (1) 

where  is  a scalar  such  that 

E * -jwA  - V$ 

For  a radially-directed  unit  current  element, 


J - J f - fS(r-r') 


. f 6(r-r'  ) 6 (6-0  ' ) g (0^4> 1 ) 


r'2  sin2© ' 


Expanding  out  (1)  yields 


-1  9 . „ ^jr  1 3 Ar  2 p 

2 90  8in  0 90  " 2.2  "775  k Ar  f 

r sin  0 r sin  0 9<p 


2 2 
9 A f , , 

, |1 r g 1 r . 

Ir  9r90l  Ir  sin  0 3 3 r ™ 


- yJrf  - jooye 


9$  £ . 1 9«t> 

9^  r s ilt  0 &4> 
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from  which  it  is 

seen  that 

3 2 A 

r 

3r39 

" -jtoue 

3$ 

36 

32A 

r 

3<$3  r 

- -jOOUE 

3 4> 
3<J> 

(4) 

These  conditions  are  automatically  satisfied  by  the  gauge 

choice  3A 

4>  - - ’rrr- (5) 


jwye 


3r 


Substituting  (2)  and  (5)  into  (3)  leaves  the  one  scalar 
component  equality 


3 2A 


2 2 

3r  r sin  0 


3 

36 


sin  0 


3 A 

i 

30 


2 

1 3 A 

r2sin  0 3$2 


+ k2A  - -y6  (r-r ' ) 
r 

which  can  be  rewritten  in  the  more  convenient  form, 

(72+  k2)(  J 

r r 

To  obtain  (7),  one  notes  that 

6 ( r - r 1 ) _ 6 (r 

r 

A solution  of  (7)  which  satisfies  the  radiation  condition  for 
a exp(ju)t)  time  convention  may  be  written  by  inspection  of  (7) 
as 


(6) 


(7) 
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(8) 
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Ur 

4ir  r ' 
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and  the  general  solution  to  (3)  for  a distributed  set  of 
currents  is  obtained  by  superposition: 


This  form  of  the  vector  potential  has  also  been  used  by 
others  [3].  For  the  symme t rica 1 ly -exc it ed  cone  and  its 
image, 

jr(r.)  . _§£- — |6(e’-e0)  - 6(e’-ir+e0)J  (io) 


where  J is  the  blcone  surface  current  density,  Sub- 
■ r 

stltuting  (10)  into  (9)  gives  the  desired  equation  for  the 
vector  potential: 


2 Tr 


A " 

r 


Ur  sin  0( 
4tt 


// 


J (r  ' ) 

sr 


-JkR 


-jkR 


dr ' 04  ' 

(ID 


where 

R~*  ^ x ^ + r ' ^ - 2 rr'  [sin  9 sin  9qCos  (4>~4> ' )+  cos  8 cos  6q] 


The  plus  superscript  denotes  source  points  on  the  upper  bicone 
surface  while  the  negative  sign  denotes  source  points  on  the 
image  surface. 


It  i 8 convenient  to  introduce  the  total  axial  current 


I (r  ’ ) * 2 tt r ' sin  0n  J (r  ’ ) 

0 s r 


(12) 


so  that  (11)  becomes 


27T  L 


■■■-Sr  I / 


e-jkR'  e-jkR 


dr  ' d<(> ' (13) 


0 0 


The  radial  component  of  electric  field  is  now  given  by 


E 

r 


1__ 

j wue 


(14) 


The  simplicity  of  (13)  and  (14)  compared  to  the  usual 
vector  potential  representations  should  be  emphasized  at 
this  point.  One  notes  that  in  the  usual  representation,  two 
vector  potential  components,  and  Aq , would  be  present. 
Furthermore,  the  integrands  of  the  potential  integrals  would 
contain  somewhat  complicated  dependences  on  angles  between 
observation  and  source  points  which  arise  from  projecting 
the  source  vector  onto  the  potential  component  vector  for 
each  source  and  observation  point,  Finally,  the  expression 
for  the  r ad i a 1 ly- d i r e c t ed  electric  field  would  be  complicated 
and  difficult  to  handle  numerically  compared  to  the  approach 
to  be  followed  here.  These  complications  Indeed  will  appear 
in  the  formulation  which  includes  a topcap  on  the  bicone  struc 
tur  e (Append ix  B) . 
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An  integral  equation  for  the  rurrent  is  obtained  by 
applying  the  boundary  condition  that  the  radial  electric 
field  must  equal  the  Impedance  loading  times  the  total 
current  density,  i.e., 


E - Z (r)l(r) 
r s 


Since  all  currents  and  fields  are  ^-independent,  it  suffices 
to  take  all  observation  points  along  the  intersection  of  the 
plane  4>“0  and  the  conical  surface.  Hence,  we  obtain  finally 


jiaye 


-C  + k2 

dr2 


A - Z (r)I(r)  - 0,  0 <r  < Jt , 6-0n,4>«O 

r s — ■ ’ 0 


(15) 


Equation  (15)  is  an  integro-dif f erential  equation  for  the 
induced  current  on  the  bicone.  As  it  stands,  (15)  does  not 
appear  to  contain  a driving  term  due  to  the  applied 
voltage  at  the  bicone  terminals.  In  the  next  section, 
however,  this  term  appears  as  a "boundary"  condition  on 


dAf/dr  at  r“0.  One  also  notes  in  (1.5)  that  discrete  or 


lumped  loading  may  be  introduced  by  allowing  Z (r)  to  be 


represented  by  appropriate  6-functions, 

N. 


V> 


L 

■ V Z,  6 ( r- r ) 
n-1  L"  L" 


t h 


for  loads  where  Z^n  is  the  impedance  of  the  n “ load 


located  at  r -r 


Ln 
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SECTION  II L 


APPLICATION  OF  METHOD  OF  MOMENTS  TO  A 
BICONE  WITHOUT  ENDCAPS 


The  usual  procedure  in  applying  the  method  of  moments 
[4]  is  to  first  represent  the  unknown  current  as  a linear 
combination  of  an  appropriate  set  of  basis  functions  and 
then  "test"  the  resulting  integral  equation  with  a series 
of  testing  functions.  Here  it  is  convenient  to  reverse 
this  order  and  to  first  test  the  equation  before  expanding 
the  current.  A set  of  testing  functions  which  offer  a 
number  of  advantages  in  a numerical  procedure  are  the  piece- 
wise  sinusoidal  testing  functions: 


w1(r) 


sin  k(Ar-r) 
sin  kAr 


, 0 < r < Ar 

, Ar  , _<  r <_  L 


w ( r ) - 

m 


sin  k f Ar- | r-r  | ) 
sin  kAr 


’ rm-l-r-rm+l 


, |r-rj  > Ar 
m-  2,3,  . . ,M  (16) 


where  Ar  ■ L/M,  r - (m-])Ar,  tn*l,2,,.,M 

m 

These  testing  functions  are  shown  in  Figure  2,  An  inner 
product  is  next  defined  as 


< f(r),  g ( r ) > - I f (r)g(r)di 


and  (15)  is  successively  tested  with  each  of  the  w , 

m 


m- 1 , 2 , . « , ,M : 


+ k2]  “?> 


Z (r ) I (r  ) , 


“■) 


0 , m-1 , 2 , . , . ,M 


Taking,  for  the  moment,  m_l  and  integrating  the  first  term 
by  parts  twice  results  in 


— + k2  A , w 
. _ 2 r’  1 


17  \(r2> 


k cos  kAf  . , . 

sin  kAr  t(T2) 


where  A (r  ) * A , Note  that  although 

r m r r " r 

ra 


r 0-TT/2 


Ar(r^)  is  not  zero  along  the  bicone.  In  fact,  one  notes  that 


'0  juiper  9r90 


and  that  the  bicone  voltage  at  r"0  is  just 


11 


E0rd6 


"<  <■  2 

1 f 3 Ar 

We  J ^10  d°  r,0 


__1_  8Ar(rl) 

jwue  Vr 


where,  using  (19),  one  sees  that  3Ar/3r“0  at  0=it/2, 
Thus  for  ro-1,  (18)  becomes 

joins  TTTTTr  [-  cos  k4r  Vrl>  + Sr(r2>l 

- <^8<r>I(r),v^>.  .V() 


For  m“  2,3,4,. ...M,  integration  by  parts  twice  in  (18)  re- 
sults in 


jU)UC  sin  kAr 


(r  t1)  - 2 cos  kArA  (r  ) + A 

UH-1  r m 


Zg(r)I(r)  , 


r<rn,-l>] 

' “a)’  ° 


m * 2 , 3 , . , , ,M 


Note  that  the  choice  of  testing  functions  has  resulted 
in  removing  all  the  derivative  operations  from  the  operator 
equations.  This  is  the  principal  advantage  of  the  testing 
functions  chosen. 


A matrix  equation  now  results  if  the  current  is 


expanded  in  an  appropriate  set  of  baei9  functions,  A 
convenient  set  is  the  pulse  functions  defined  by 


Pj(r>  - 


Pn(r) 


1 , 0 < r <_  Ar/2 

0 , Ar/2  < r < L 


1 , | r-r  | < Ar/2 

n — 

0 , | r-r  | > Ar /2 

i n i 


n - 2 , 3 , , , , M 


(see  Figure  3)  and  the  resulting  current  expansion  is 


M 


I(r)  3 l.  Inpn(r) 
n-1 


(22) 


Note  that  the  current  at  the  bicone  edge  reL  is  auto- 
matically zero  by  our  choice  of  basis  functions  (Figure  3), 
When  (22)  is  substituted  into  (20)  and  (21),  there 
results  the  system  of  linear  equations 


June 


“sitTkAr  ['  cos  kAr  4-  (ri;ri,r1  + ) + ¥ ( r 2 ; r 2 , r l + >] 

- ^za(r)Pi(r>  * wi 
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' 1 
; SJ 

; 

i ? 


. — - -M/J4.U*. v wl. it tflfcaaAtjy 


\w*\> 


y x < — ■ — — 

L n ' j U)\i  r sin  k 


L'  CO‘  kA''  nrl;ri,--rnt!  + 'l'<r2:r«-'rn  + ). 


-<Z* 


( r)p  (r),  i»  1 >>-  -V 


0 

(23) 


and 


1 ] jojue  sin  kAr 


jj*,('rm-l;rl,rl  + )_  2 COS  kAr  y(rm;  rl  ,r2^ 


+ Y(r 


m+1 


:rlr3  + 0 


‘ I I \ 7—: — . 7-.c  [V(r  , ; r _ , r + )-2  cos  kAr'V(r  ;r  _,r  +) 

n.T  n \ jojyc  sin  kAr[_  m-1’  n * nT  m’  n n+ 


+ ^(rm+l;rn-*rn  + >]  " \Zs(r)V(r)’W 

» 0 ,ra«2 , 3 , . . ,M 
(24) 


These  equations  may  be  assembled  into  the  matrix 


equa  t ion 


where 


ZI  - V 


\ 1 

“-Vo" 

0 

1 = 

, V - 

0 

0 
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and  the  elements  of  the  impedance  matrix  may  identified 


f rom  (23)  and 
defined  by 


¥ (r  ; r - , r 4.) 

n nT 


(24).  The  functions  V(r;r  ,r  + ) 

n-  nT 


are 


d r ' d 4> ' 


(25) 


where 


R’ ■ r^  + r 1 ^ - 2rr'(8in^GQ  cos  <j) ' +cos^Bg)  (26) 

and  where 

r ± - r + Ar/2 
n n - 

The  evaluation  of  the  double  integral  in  (24)  is 

simplified  by  analytically  approximating  the  integration 

with  respect  to  r'.  This  may  be  accomplished  by  noting  that 
in  general  k(r  +-  r ) <<  1,  so  that  a few  terms  <•  a 
Taylor  series  expansion  about  some  point  rn  in  the  interval 
[rn_,rn4_j  should  be  accurate.  Ac  c o r d i n g 1 y,  one  writes 

-JkR*  - j kR4  -jk(R4-R4) 
e * e e 

- 1 k R ~ + + 

- e " i 1 - i k < R “ - R~)]  (27) 

where  R~“  ^~lr»r  * Substituting  (26)  into  (24)  and  listing 
n 

as  a fourth  argument  the  point  about  which  the  expansion  is 
made,  one  obtains 

1 6 


I 


* -r--  ••  , ...»  ... 


V(r;r  ,r  ) ~ f(r;r  ,r  r ) 
n-  n+  n - n+  n 

2 n r + 


8ti 


r f {"  [-><  Ik 

2 J J L t''r’R+  ' 


~ j kR, 


-e 


1+jkR 


B 1L 

r ' R"  r ' 


. JiL 

4 TT  2 


* 

r 

/ i 

0 

-J“n 

e 

L 

i+jkR; 

n , 

r log 

rn+<C  + r ' rn-b+) 


r (R+,  + r - r b+) 
n-  n+  n+ 


dr  ’ d<t 


- j k£n 


r + 
n 


r _ 
n 


* J kR 


-e 


1 + jkR" 


log 


• +(R  _ -t-  r - r b ) 

n+  n n 


rn-(Rn+  + r ' rn+b’> 


-j  k«,n 


rn  + 

r 


n _J 


d<J>  ' 


(28) 


+ 2 ,2 

where  b ■ sin  0^  cos  $ + cos  0^.  In  several  situations. 


appropriate  limits  of  the  integrand  of  (27)  need  to  be  taken 
First,  when  the  source  is  the  current  segment  at  the  bicone 
terminals,  the  integral  (27)  reduces  to  the  simple  form 


•Mr  ; r j , r 1+) 


^*  ( r 5 r 2 » r j + > r j ^ 


4 TT  ‘ 


f -Jkr 

1 e (l+)kr)log 

R”+  + r - b"r1+ 

+ . -f 

/ 

0 

R,+  + r - b r 2 + 

d <J> ' 


(29) 


1 7 


'i.U. 


As  the  observation  point  r in  (28)  approaches  the 


bicone  terminals,  r*-r^*  0,  the  limiting  form  of  the  inte« 
grand  can  be  integrated.  The  result  is 


4'(rl;rl,rl+,rl)  “ 2tT  lo8(tot  T5 


* 

□ L 

f I 


a very  interesting  result  that  is  independent  of  the  sub- 

domain  size  at  the  bicone  terminals.  Finally,  all  the  so- 

called  "self  terms"  Y(r  ;r  ,r  x , r ).  ni*  1 , contain  an 

n n - n"r  n 

integrable  singularity,  In  fact  one  easily  establishes 


r (R+  + r - b+r  _) 

log  — “ ~ 2 lo8  U ' I 

r _ ( R , + r - b r +)  rn 
n n+  n+ 


This  singular  term  is  then  subtracted  from  the  integrand  in 
(27),  resulting  in  a non-singular  Integrand  which  is  then 
numerically  integrated.  The  term 


H n | <J) 1 ] d o ' “ y^-  ( i,  n it  - 1 ) 


is  then  added  to  the  result  to  take  care  of  the  part  of  the 
integral  contributed  by  the  singularity, 


FORMULATION  AND  NUMERICAL  SOLUTION  OF  AN 
INTEGRAL  EQUATION  FOR  A BICONE  WITH  ENDCAPS 


The  formulation  c f the  integral  equation  for  a cone 
radiator  over  a ground  plane  with  an  end  cap  is  considerably 
more  complicated  than  that  for  the  case  when  the  endcap  is 
not  present.  It  Is  possible,  however,  to  generalize  the 
approach  used  for  the  bicone  without  a topcap  and  to  trans- 
form the  derivatives  appearing  in  the  equations  into  har- 
monic operators  along  the  radial  cone  and  topcap  coordinates, 
as  is  done  in  Appendix  B.  This  approach  has  the  advantage 
again  that  testing  with  piecewise  sinusoids  allows  the  re- 
placement of  derivatives  by  a finite  difference  of  potentials. 
However,  to  effect  this  transformation,  an  extremely  com- 
plicated kernel  must  be  used  (see  Appendix  B ) which  contains 
many  singularities  other  than  the  usual  ones  where  source  and 
field  points  coincide.  While  this  approach  has  been  used, 
it  has  been  found  to  be  unwieldy  and  rather  inefficient. 

The  approach  described  here  begins  with  the  description 
of  fields  in  terms  of  the  more  commonly  used  vector  magnetic 
and  scalar  potentials  expressed  in  terns  of  the  cone  currents 
and  charge.  However,  it  is  Found  that  these  potentials  are 
singular  at  the  bicone  terminals  which  again  creates  an 
unnecessary  complication.  In  order  to  circumvent  this  problem, 
the  cone  and  Image  surfaces  are  allowed  to  Intersect  with  a 


h’r’VM.r- 


small  "waist"  of  radius  "a"  (Figure  4).  If  "a"  is  very 
small,  there  should  be  negligible  difference  in  the  input 
impedance  and  currents  found  for  this  case  and  that  for  the 
limiting  case  of  a *0 . For  convenience,  the  cone  coordinates 
are  defined  with  respect  to  the  projection  of  the  cone  sur- 
face to  a tip,  as  in  Figure  h.  Furthermore,  the  direction 
of  the  unit  vector  r and  the  positive  direction  of  corres- 
ponding vector  components  is  taken  to  be  towards  the  center 
of  the  topcap,  in  the  direction  of  decreasing  rt* 

The  integral  equations  are  obtained  by  setting  the 
radiated  field  tangent  to  the  cone  surface  equal  to  the 

impedance  drop  per  unit  length  due  to  the  loading: 
s <6 

-jwAr  - - ZsIc  - 0,  a/sin  0Q<  rc<  L + a/sin  0Q  (31) 

c c 


- J <*)A  - - Z I - 0,  0<  r < L sin  0_  + a (32) 

r ^ **  c sc  t — u 

where  the  tangential  components  of  magnetic  vector  potential, 

A and  A , are  given  by 
rc  r t 


2 tt  L+a / sin  0 


A “ 

p 8TT 


or  F + -JkRf 

c 08  t,  e 


0 a/'sin  9 , 


m — ^ 

v pc 


-jkR 

cos  C “e  PC1 

+ ^ dr’< 

R I c 

pc 


2tt  Lsin9_+a  , -jkR 

» » Or  t + pt 

f r cos 


II 

0 0 


, -jkR  „ 
c + pt' 

cos  f,  e 


dr  ^ d 4>  ' 
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The  scalar  potential  is  given  by 


I,+a  / sin  8 , 


V / / 


a / s l n 6 , 


- i k R - i kR 

e Pc  e ‘ pc 


dr  d<t> ' 
c 


2tt  k sin  0 n + a 

Am  ” 


/ / 

0 0 


d I 

t 

(.'K. 

-jkR 

e 

d r t ‘ 

R\ 

Pt 

dr^dij)  ’ 


The  currents  1 and  I are  the  total  linear  currents  on 
c t 

the  conical  and  topcap  surfaces,  respectively,  and  are  re- 


lated to  the  corresponding  surface  current  densities  and 


J by 
t 


I =2ttr  sin  0nJ 
c c 0 c 


I * 2nr  J 
t t t 


The  distance  quantities  are  all  of  the  form 


+ / ' 2 + ' + 
R ~ « r + 2 b ' r +c 

pq  q pq  q pq 


p , q - c , t 


where 

+ 2 2 _ 

b » - r sin  0.  cos  $ 1 + r cos  6 - a cot  9. cos  0_(1+1) 

cc  c 0 c 0 00 
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1 vrja^'-^iT'- 


, 2 2—2  — . 
r = r ‘ + a (’of  <\(1  + 1)  - 2r  a cos  B cot  <i  (1-t-L) 


c i c 


0 0 


b “ «-r  sin  9-  cos  A* 
ct  c 0 


+ 2 2 
c - r - 2r  cos  0A(a  coc  8.+  L cos  6-)  + (a  cot  0n+  L cos  0) 
etc  c u u—  0 U—  u 


b " «+L  cosl„  - a cos  9Acot  BA  - r sin  0A  cos  4> f 
tc  0 OOt  0 


c ' ■ r^+(L  cos  0A  + a cot  0 
t c t U — 0 


btt  “ -rt  cos  <|>  * 


+ 9 o o - 2 

c 1 1 - r t + L cos  0o(l+1) 


The  angles  between  the  source  current  elements  and  the  tan- 
gential component  of  electric  field  at  the  observation  point 
are  determined  by 


+ 2 2 
cos  £ - + cos  d>  ’ sin  9„  + cos  8_ 

’cc  0 0 


cos  £ 7 m + sin  0,,  cos  $ ’ 
ct  0 


CO  8 £ “ * + COS  * 


cos  £ ■ + cos  sin  0. 


It  is  convenient  to  choose  as  testing  functions  the 
pulse  functions  shown  in  Figure  5.  Thus,  testing  (31)  with 
p ^ results  in 


<Ar  *P1>-  <jv-’  Pi)  " <Vc  *pl)  “ ° 


Upon  integrating  by  parts  in  the  central  term,  and  noting 

that  A is  slowly  varying  over  the  interval  and  hence  may 
c 

be  approximated  by  (r  ^),  one  obtains 

c 


-J“\  '■'cl’  — ‘ l(4<rcl+  4ro/2)  - *<rcl 

c 


>1  - <Vc  -Pi)-  0 


But  4>(r  j)  is  just  the  bicone  terminal  voltage  VQ  with  res- 
pect to  the  ground  plane.  Hence, 


-jujA  (r  .)Ar  - 2 4>  ( r .+  ir  12)  - /z  I , p.\ 
rc  cl  c cl  c \ s c 1/ 


For  the  remaining  testing  functions  on  the  cone, 
testing  of  (31),  integration  by  parts  on  t-he  scalar  potential 
term  and  approximation  of  the  vector  potential  by  its  value 
at  the  center  cf  the  pulse  yields 


■jwAr  A (r  ) - [ 4>  < r + Ar  /2)  - ( r - Ar  / 2 ) J 

c r cm  cm  c cm  c 


- <Vc  . P.>  - 0.  ■ - 
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At  the  edge,  the  testing  pulse  consists  of  two  parts, 
one  on  the  cone  surface  and  one  on  the  topcap.  Hence,  both 
equations  (31'  and  (32)  must  be  used.  With  approximations  on 
the  two  vector  potential  componants  similar  to  that  above, 
integration  by  parts,  and  enforcement  of  continuity  of 
the  scalar  potential  at  the  edge,  one  obtains 

rArc  Art  1 

hr  Ar  (rcN  > + ~r  Are(rtl)J 

-t»(r  -it  /2)-«<r  +4r  /2U  - P,  > - 0 (36) 

c ' c ' 

where  1^  ■ Ic  or  I ^ as  is  appropriate.  On  the  topcap, 

one  has,  analogous  to  (37), 


■J“4rc*rt<rc.)  - itt/2)  - *<rtm+  it£/2)] 


<Vt  • PN  -1*„>  ' °-  ""2 V1  (3,) 


The  current  is  next  expanded  in  the  set  of  pulse  functions 


Pn  of  Figure  3, 


N 


'c  <rc>  * l 

n>*  1 


(40) 


N + N 
c t 


K (r,.>  “ l IP  (r.) 
t t nnt 


(4  1) 


Note  that  the  current  at  the  edge  of  the  cone  is  the 

c 

same  on  both  the  cone  and  the  topcap  surfaces.  The  derivatives 


of  the  currents  above  are  approximated  by  a finite  difference 
of  adjacent  current  pulses  which  is  then  assumed  to  be  ex- 


panded in  its  own  set  of  "charge"  pulses  (see  Figure  6 ); 


dl  Nc-1  fl  - I 

c ..  V n+1  n + , , 
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aT  N fN  (.  _ 
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177  * ”477'~  p»<rt> 

t n-N  t 

c 


where  1^  +1  is  taken  to  be  zero, 

c t 


Thus  the  vector  potential  quantities  in  (3f>)-(39>  may 


be  written  as 


i - -S  (l, 
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The  inner  integrals  In  (46)  and  (4  7)  ran  be  approximately 
analytically  integrated.  They  are  all  of  the  form 


r 

/ 


n + -jkR 

e 


d r ’ 


(4  8) 


n - 


where  R is  of  the  form 


vr*2+  2r'b  4 c 


(4  9) 


Since  the  range  of  integration  [rn_»  rn_jJ  8maH  com- 


pared to  a wavelength,  it  is  appropriate  to  expand  e 


■jkR 


in 


a Taylor  series  about  some  point  which  is  the  distance 
from  the  observation  point  to  a point  rn  in  the  interval 

frn-»  rn+]-  Thus’ 


-jkR  -jk(R-Rn)  -jkRn 
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The  error  in  the  real  and  imaginary  parts  is  less  than 


Max 
r ' e [ r - , r 
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n ’ n + 
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where  Ar  is  the  subdomain  size,  For  five  subdomains  per 
wavelength  (Ar/X  ■ 1/5),  this  results  in  a maximum  error 
of  less  than  IX  in  both  the  real  and  imaginary  parts  of 
the  integrand.  The  resulting  integral  should  indeed  be 
much  more  accurate  than  this.  With  this  approximation, 


n+  -J kR 


/ * 


- j kR 

dr'  * e n 


(Ir  JI2) 
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Combining  these  results,  one  has,  finally, 
-fn+  -JkR 
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Equation  (50)  is  singular  in  <p ' if  the  observation  point  is 
in  the  interval  [rn_,rn+l.  Hence,  the  integrals  in  (4  6)  and 
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(47)  need  to  be  evaluated  by  subtracting  the  singularity 

from  the  integrand  and  adding  its  integral  to  the  numerically 

determined  integral.  If  r is  in  the  interior  of  the  inter- 

n 

val,  r <r  < r , (50)  behaves  like  — 2 Sin  I cb  ’ I near  d> ' - 0; 
n-  n n+  11 

if  r ■ r or  r - r , (5  0)  behaves  like  - Sin  I <J> f I . The  de- 
n n-  n n+  1 ' 

tails  of  the  procedure  parallels  that  described  at  the  end 
of  Section  III. 
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SECTION  V 


NUMERICAL  RESULTS  AND  CONCLUSIONS 


This  section  describes  numerical  results  obtained 
from  the  computer  code  developed  from  the  theory  described 
in  Section  IV.  The  resulting  code  was  written  to  model  the 
cone  with  or  without  a topcap.  Hence,  results  from  this 
general  code  could  be  checked  against  those  obtained  from 
the  code  based  on  the  methods  of  Section  III  for  the  cone 
without  a topcap.  For  narrow  cone  angles,  the  calculated 
input  impedance  for  unloaded  cones  for  various  frequencies 
was  also  compared  to  the  theory  of  Schelkunoff  [5l  and 
found  to  be  in  very  good  agreement.  Results  from  the  general 
code  foi  moderate  cone  angles  were  also  compared  with  those 
computed  by  the  method  of  Appendix  B,  which  includes  the 
effects  of  the  topcap.  These  comparisons  were  made  to  vali- 
date the  consistency  of  the  various  approaches  and  to  com- 


pare with  existing  data.  It  was  also  established  that 
the  input  reactance  at  low  frequencies  could  be  used  to 


check  the  atatic  capacitance  calculated  in  the  companion 
report  [1]  for  both  the  loaded  and  unloaded  case.  Finally, 


it  was  verified  that  the  computed  results  were  almost 
independent  of  the  choice  of  the  waist  radius,  a,  of  Section 
III,  provided  a was  chosen  small  enough. 


jki  .L 


All  the  data  In  this  section  pertain  to  a conical 
antenna  with  a vertical  height  of  40  meters  and  a cone 
angle  of  G^**  42.26°.  Theae  parameters  translate  to  a cone 
slant  height  of  54.05  meters  and  correspond  approximately 
to  the  cone  considered  in  [6l.  The  locations  and  valuee 
of  the  lumped  resistive  loads  used  are  listed  in  Table  1 
and  are  taken  from  [6]. 

Figures  7-10  illustrate  the  current  distribution  on 
the  cone  at  a frequency  of  825  KHz,  approximately  the  first 
resonant  frequency  of  the  unloaded  structure.  Figures  11-14 
illustrate  the  same  results  at  1.375  MHz,  approximately 
half-way  between  first  and  second  resonance  of  the  unloaded 
structure  (see  Figures  15  and  16).  Two  features  of  the  cur- 
rent distributions  are  notable.  First,  the  edge  condition 
[7],  which  requires  that  the  current  at  the  edge  has  in- 
finite slope,  and  the  continuity  equation  relating  current 
and  charge,  which  requires  that  the  total  current  approach 
zero  with  zero  slope  at  the  center  of  the  topcap,  combine 
to  limit  the  amount  of  current  the  topcap  can  support. 
Secondly,  the  loading,  which  increases  to  a maximum  at  the 
edge,  further  limits  current  flow  on  the  topcap. 

Figures  15-18  illustrate  the  variation  with  frequency 
of  the  input  impedance  of  the  conical  structure  for  the  vari- 
ous loading  and  topcap  configurations.  Again,  the  influence 
of  the  topcap  is  found  to  be  negligible.  The  absence  of 
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Table  1.  Positions  and  values  of  loading 
resistors  on  the  cone* 


i 

5 

< 

V 


LENGTH  ALONG 
CONE  GENERATOR 

RESISTANCE 

(METERS) 

(OHMS) 

12.17 

4 .69 

14.34 

7.17 

16.90 

9.06 

19.83 

11.74 

23.54 

15.43 

27.73 

20.82 

32.59 

30.36 

38.41 

50.38 

45.31 

114.88 

53.43 

114.88 

63.15 

100.00 

72.25 

100.00 

81.35 

100.00 

CURRENT  *10' 


f (MHz) 


Figure  18.  Input  Impedance  of  loaded  cone  without 

topcap,  L“54.05m.  0q«  42.26°.  Imaginary 
values  are  negative. 
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resonances  .In  the  loaded  case  can  be  attributed  to  the 
effectiveness  of  the  resistive  loading  in  eliminating 
reflections  from  the  cone  edge  which  would  result  in 
standing  waves  on  the  structure. 

Radiation  patterns  in  the  near-field  region  of  the 
loaded  structure  with  a topcap  (r  - 100  meters  and  the 
frequency  is  550  KHz)  are  shown  in  Figures  19-21.  Figures 
22-24  give  the  corresponding  patterns  in  the  far  field 
(r  ■ 10^  meters).  For  comparison,  far  field  patterns  for 
the  unloaded  structure  are  illustrated  in  Figures  25-27. 

Although  resources  did  not  permit  a time-domain  analysis 
of  the  response  of  the  structure,  auch  a study,  which  could 
include  a simple  equivalent  circuit  model  of  the  puiser, 
would  be  a logical  extension  of  the  present  problem.  To 
be  done  efficiently,  however,  some  improvements  In  the 
present  computer  code  should  be  implemented.  Specifically, 
an  adaptive  Integration  procedure  should  be  employed  to 
handle  the  integrations  over  the  conical  current  sub- 
domains,  whose  radii  vary  drastically  from  regions  near  the 
feed  to  those  near  the  cone  edge.  The  present  code  uses  a 
fixed  order  quadrature  rule  for  all  segments  on  the 
structure.  Additional  parameter  studies  can  be  carried  out 
using  the  present  code  to  assess  the  effects  of  lumped  vs. 
distributed  loading  and  the  effects  of  various  load  distri- 
butions on  the  performance  of  the  simulator.  A more 
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Figure  19.  Efi  radiation  pattern  for  loaded  cone 
with  topcap,  L-54.Q5m,  0O*  42.26°, 
V_»l  Volt,  f - 550  KHz.  r - 100m. 


? 

f 


Figure  20.  H.  radiation  pattern  for  loaded  cone 
with  topcap,  L*54.05m,  9 “ 42.26°, 
V0«  1 Volt,  f - 550  KHz,  r - 100m. 
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Figure  23.  H 


w^th  topcip"  io*ded 

V - 1 Vn,+  « .:U3°*  9n"  42.26®. 


V - i u i " J'4  • u-,o , 

v0  1 Volt,  f - 550  KHz 


90“  42.26'  , 

r " 10*m! 
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Figure  27.  E radiation  pattern  for  unloaded  cone 
with  topcap,  L«*54.05m,  0 - 42. 26°, 

V - 1 Volt,  f - 550  KHz.  r - 104m. 


ambitious  project  would  more  carefully  model  the  actual 
wire  structure  with  loading. 


One  conclusion  of  this  and  the  companion  study  [l] 
is  that  the  addition  of  a topcap  does  not  significantly 
change  the  electromagnetic  parameters  of  the  structure  - 
at  low  frequencies,  the  static  capacitance  and  effective 
heights  are  almost  unchanged  and  at  the  higher  frequencies 
the  loading  and  the  sharp  angle  at  the  edge  tend  to  pre- 
vent current  from  flowing  on  the  topcap.  This  observation 
may  have  some  impact  on  the  design  of  future  simulators. 
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APPENDIX  A 

CALCULATION  OF  RADIATED  FIELDS 

Once  numerical  values  far  the  current  distribution 
have  been  determined,  the  fields  radiated  by  the  bicone 
structure  can  readily  be  determined.  Because  of  the 
symmetry  of  the  structure  and  the  excitation,  the  only 
non-zero  components  of  the  electric  and  magnetic  fields 

are  E , En , and  H,.  These  are  defined  in  terms  of  the 

r t)  $ 

vector  and  scalar  potentials  as 

Ef  . -i„At  - 


Ee  - -J»A6  - i f| 

. 3Afl  - 3A 

uH.  ■ —A „ ■■  — w ijl > — L. 

H 4>  r e 3r  r 30 

A A 

where  A - A^r  + A^B.  A spherical  coordinate  system  centered 
at  the  b leone  feed  and  with  9 measured  from  the  z-axis  is 
assumed.  Since  the  fields  are  (((-independent  , all  fields 
are  evaluated  in  the  x-z  plane  where  <fs  — 0 . The  vector 
potential  and  scalar  potential  are  given  by 
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8tt  jure  n-1 


Ar  rc 
c ' 


c 1 -I  Ar 

I -4~  t.e-V  nr 


N I ,,-I  Ar 

+ I +1  -Trr  ♦«  ~r 

n-N  +1  t 

c 


(A-2) 


(A-3) 


where  the  currents  and  coordinates  are  defined  in  Section  IV. 
The  potential  functions  ¥ and  are  defined  as 


Y (r,0,r  )-Ar 

pq  q q 


2 IT 

/h"*? 

f I -><» 

■"'/  h 


+ cos  £ 


d*'  (A- A) 


r i -K  -Jk% 

<P <r,9,r  )-Ar  j [ - ■£ 

q q q J 1 Rt  R" 

0 \ q q 


d(j)' , p-r,0;  q-c  , t (A- 5) 


where 

cos  fj  * • ± sin  8 cos  <t>'  sin  8„  + cos  8 cos  8. 
rc  0 0 

cos  ■ i cos  8 cos  *'  sin  8_  - sin  9 cos  8„ 

8c  0 0 

cos  * + sin  9 cos  <J>' 

cos  £g*  * + cos  9 cos  ij)' 

and  the  radius  vectors  are  all  of  the  form 


R / 1 2 . + , . ± 

q r + 2b  r + c 

q q q q 
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f 


with 

t _ 

» -r  cos  sin  9 sin  00  +■  r cos  9 cos  00 

±2  2 2 
c ■ r ± 2ra  cos  6 esc  9 + a esc  0„ 
c 0 0 

+ 


- a cot  9, 


■r  sin  0 cos  <j> ' 


c*  - + 2rL  cos  9 cos  6.  + 

c 0 


Assuming  a suitable  choice  for  Ar  and  A9,  one  may  approxi- 
mately compute  the  fields  in  (A-l)  by  finite  difference 
approximations ; 


Er(r,0)  - -juAr 


(r  6)  - *lr±Lr.,  6) 

(t’0)  Ar 


Eg (r , 9 ) 


“jwAg (r , 9) 


-H laltAg ,6) 

rA0 


^A0(r,9)  + 


AQ(r+Ar ,9)-Ag(r ,9) 
— 


A^ (r ,0+A6)-Ar (r  ,9) 

Tab 


(A-6  ) 
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APPENDIX  B 


AN  ALTERNATE  INTEGRAL  EQUATION  FOR 
A CONE  WITH  TOPCAP 


The  purpose  of  this  appendix  is  to  show  how  a novel 
identity  involving  the  free  space  Green's  function  may  be 
used  to  change  the  integral  equation  into  a form  where  the 
testing  procedure  of  Section  III  is  applicable.  In  ex- 
change for  simplicity  in  the  form  resulting  from  the  testing 
procedure,  however,  one  obtains  extremely  complicated  ker- 
nels in  the  integral  equation.  Furthermore,  the  new  ker- 
nels have  a number  of  singularities  other  than  the  usual 
one  where  match  points  and  field  points  coincide.  These 
complications  make  both  the  analysis  and  the  numerical 
treatment  tedious.  Nevertheless,  numerical  results  have 
been  obtained  for  several  cases  using  the  approach  and  the 
results  are  iti  good  agreement  with  data  obtained  by  the 
method  of  Section  III.  For  simplicity,  we  treat  here  only 
the  unloaded  cone. 

As  a prelude  to  the  integral  equation  derivation, 
we  derive  a transformation  of  the  formula  for  electric 
field  components.  Consider  the  x-component  of  electric 


field  given  by 


jcopE:  Ex  * x*(k  +VV»)A 

“ k2v 


2 2 

~2  9 3 A 3 A 

( — + k ) A + ^ — 5-^-  + ~ — ~~ 
-2  x 3x3y  3x3z 

d X 


( B-  1 ) 


where  the  vector  potential  A in  terms  of  current  density 


J is 


A » 


Ait 


/ 


~5k\r-r'\ 


j 5_ 


dV’ 


(B-2) 


The  identity 


3 u3  v 


-jkR 


J_  + k2 

3u2 


uv 


fe-jkR 


2 2 
v 4-w 


(B-3) 


/ '2  2 2 

where  R - u +v  +w  can  be  used  with  u ■ x-x',  v - y-y'. 


and  w ■ z-z'  to  rewrite  (B-l)  as 


j (x~x  ' )_(y-y  ' ) 
y (y-y ' )2+(x-z  ' )Z 

-(x-x'Hgl?'? dV' 

(y-y ' ) 2+(z-z ' ) 2 J R 


) [9  (y-y  ' )-tz  (z-z ' ) ] 
y’)2+  (z-z’)2 

-jkR 

x 2-j — dV'  ( B-4 ) 

The  vector  r-r'-  (x-x')x  + (y-y')y  + (z-z')z  can  be  written 
as  the  sum 

r-r'  - (r-r')x^  + (r-r')tf 

where 

(r-r')^  " 5c*(r-r’)  - x-x* 
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I 


is  just  the  component  of  r-r'  along  the  direction  of  x and 


( r-r  ' ) ^ 6 - r-r'  - (r-r ' ) ^8 

is  just  the  component  of  r-r'  transverse  to  x.  Thus, 


(B-4)  can  be  written  as 


June  E - + 

X ^ * 


_ / £(r-r')  (r-r')  \ -jkR 

j.  s -tV-—--  ~~  dV’ 

l (r-r')J 


( B-  5 ) 


Since  the  choice  of  the  coordinate  system  is  arbitrary,  we 
may  choose  the  x-axis  parallel  to  some  constant  unit  vector 
a and  write  the  component  of  electric  field  in  the  direction 


of  a to  be 


jwye  E, 


■(b  *>■]/'•(• — ,tttH 


l/’i' 


- - le'J 
C cot  (a, r-r')  


-jkR 

HV  ' 


(B-  6 ) 


where  now  £ denotes  the  direction  of  the  component  of 


(r-r')  transverse  to  a and  s denotes  distance  along  a line 
i the  direction  of  a.  Note  that  the  integrand  in  (B-6)  is 
singular  not  only  when  R"0,  but  also  when  the  angle  between 
a and  r-r'  becomes  either  0°  or  180°.  Since  in  (B-6)  the 
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only  differential  operator  la  the  harmonic  operator,  then 
along  a line  in  the  direction  of  a,  the  testing  procedure 
of  Section  III  which  uses  piecewise  sinusoids  may  again  be 
used  to  transform  the  harmonic  operator  into  a finite  dif- 
ference operator.  This  Is  the  advantage,  gained  at  the 
expense  of  obtaining  a more  complicated  kernel,  of  em- 
ploying the  transformation  (B-3). 

Returning  to  the  cone  problem,  we  choose  the  direction 
of  S to  be  along  the  cone  generator  formed  by  the  inter- 
section of  the  4>*>0  plane  and  the  cone  surface,  and  apply 
the  boundary  conditions.  After  some  straightforward  but 
tedious  vector  projection  operations,  one  arrives  at  the 
Integral  equations 


1 

_ 2 -N 

» + k2 

• 

¥ + ¥ 

ju>ye 

3 r 

c 

cc  c t 

. . 

- V05(rc-0+>, 


1 

a2  2 

° + k 

¥ + ¥ 

jojye 

a 2 
13  r. 

tc  tt 

k , 

0 < r <1, 

— c — 


“ 0,  0<_r^<^L  sin  0Q 


(B-7  ) 
( B-8  ) 


where  2 tt  L 

¥ (r  ) - I I I (r  ’ ) (K+  + K"  )dr'd$' 

pc  p 8tt 2 J J c c pc  pc  c 

0 0 


V<V‘  a„2 


2 tt  Lsin0, 


II 

0 0 


VrP(KJt+ 


P * c or  t 


(B-9) 
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The  kernels  K~  in  (B-9)  are  of  the  form 

pq 


K (r  ,r',4>')- 

pq  p q’ 


+ C4  ( r ’ , <{> ' ) [ A±  (<f>’)r'+B4  ($')] 

D4  ($.)+  -EJ3 3 -L-LSL 1 a E i 


pq 


R42  - C±2(r’  ,$') 

pq  pq  q 


-jkR 


pq 


pq 


The  distance  between  source  points  and  field  points,  , 


is  of  the  form 


R4  - ^ r ' 2 + 2b4  r ' + c4 


PT 


pq  q pq 


+ + 

and  b and  c~  are  as  defined  in  Section  IV  with  "a"  set 

pq  pq 

equal  to  zero.  The  term  may  be  expressed  as 


C4  - e4  r'  + f4 
pq  pq  q pq 


The  parameters  A4  B*  , D4  , and  f+  are  defined  in 

pq  pq  pq  pq 


Tables  (B- 1 )- (B-4 ) . 

Testing  Equations  (B-7)and  (B-8)  with  piecewise  - 
sinusoidal  testing  functions  as  in  Section  III  results  in 
the  equations 
k 


{- 


jaTyealnkAr  \ C°8k4rc^cc(rcl)  + ^ct^cl^ 


+ ( r ,)  + Y „(r  ,)]  ) - -V. 

1 r.  p r 7 r t r 7 1 f] 


(B 


n 


A 


± 

cc 


B 


i 

cc 


D 


± 

cc 


f 


± 

cc 


TABLE  B-l 

DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  K 

cc 

sin2  0Q[sin2<t>’  + cos20q(cos  $’  + l)2] 

0 

±2  2 
-e  ■ sin '0A  cos  $ ' + cos  0n 
cc  0 — 0 

r 

c 


TABLE  B-2 

DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  Kct 

A * - cos20_  cos2ifr'  + sin2<t>' 
c t u 

i — 2 

B . * + L cos  sin  6n  cos 

c t U u 

+ 4 A * 

D “ ■ - e ^ - sin  0 n cos  <j> ' 
c t c t 0 

f ‘ * r + L cos2©., 
etc  0 


1 


TABLE  B-3 

DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  K 


tc 


tc 

± 


8in^0Q  sin^’  + cos^0q 


B . - - L cos  0, 


tc 


tc 

± 


- e * sin  0rt  cos  d>  * 
tc  0 T 


f - r 


tc 


TABLE  B-4 

DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  K 


tt 


tt 

± 


sin^$  * 


B ‘ - 0 


tt 

± 


tt 
. * 


®tt  “ coe 


tt 


■ r . 
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and 


{ 


t 

1 

I 


jwyesinkAr"  ^cc^c.m+l'  + t ( rc  , m+1  ^ 


2coslcAr  [V  (r  ) + ¥ (r  )] 
cLcc  cm  c t cm 


+ l¥  (r  ,)  + ¥ fr  . ) ] \ - 0, 

1 cc  e,m-l  ct  c , m-1  J f 

m-2 , 3 , . . . ,Nc-l 
(B- 1 1 ) 

Testing  at  the  cone  edge  with  a piecewise  sinusoidal  testing 
function  which  straddles  both  the  cone  and  the  topcap  and 
which  has  its  peak  value  at  the  cone  edge,  one  obtains 

Jii)ueolnk»tc  c0*  l4rtl,tc,rt»c)  + ’ct^cNc^ 


4 jo.iie.i«k4rt  {-  C0‘  k4rt!','tc(rtl)  4 ytt(rtl)l 
4 IV"t2>  4 »tt(rt2>>} 


1 , 

f 

3 v 

cc  , c t 

,L 

"®ftc  "tt 

joipe 

k. 

3r  3r 

c c 

V — 

T 

r -r  M 
c cNc 

3rt  *rt 

( B — 1 2 ) 


Finally,  testing  on  the  topcap  surface  yields 

jl*HJEslnkArt  4 y t t ( r t ,m+l ^ 

- 2 co.  Urt[Ttt(,[t)  + »ft<* 


tm 


)]  + 


1 

J 

i * 

J 

'I 


1 


*4 


-1 
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m-2,3, . . . ,Nt+l  (B-13) 


Substitution  of  the  current  expansions,  Eqs,  (40)  and  (41)  of 
Section  IV,  into  (B-10)-(B-13)yields  a matrix  equation  for 
the  determination  of  the  unknown  current  coefficients.  Be- 
cause of  the  pulse  expansion  for  the  current,  the  matrix 
elements  involve  integrals  like  (B-9)  but  with  the  current 
in  (B-9)  equal  to  unity  and  the  limits  on  the  radial  inte- 
gration replaced  by  the  limits  of  the  corresponding  current 
subdomain.  According  to  (B-12),  the  term  at  the  edge  also 
requires  the  derivative  of  such  integrals.  In  the  following, 
we  present  a procedure  for  approximately  evaluating  the 
radial  integration,  leaving  the  integration  to  be  done 
numerically.  The  required  integrals  are  all  of  the  form 


' n+ 


' n+ 


-jkR 


dr 


( B-l 4 ) 


n- 


n- 


where,  for  convenience,  all  subscripts  and  superscripts  have 
been  suppressed.  Since  the  number  of  subdomains  should  be 


cnoaen  such  that  k(r 


n+ 


rn_j  is  small,  we  chose  some  point 


r in  the  interval  [r  ,,r  ] and  expand  exp(-jkR)  in  a Taylor 

n-f  "* 


series  about  the  point  r1*  r^; 
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-jkR  -jkRn  -jk(R-Rn) 
e * e e 1 


-jkR 

“ e n [ 1- j k (R-R  ) ] 


(B-l 5 ) 


where  R^  denotes  R evaluated  at  r'-r^.  The  resulting  approxi- 


mate integral  is 


1 n+ 

/ ■ 


dr'  * e 


'•ft-  ^[^]- 


-jkR 


(I1  + 12  + I3) 


where 


- n+ 

l"3kJ  1° 


R -C 


(B-l 6 ) 


12  - DU+jkR 


■>/ 


13  - (1+jkR 


•'/ 


r2-c2 


The  first  integral  may  be  evaluated  by  substituting  the 
definitions  for  R and  C in  terms  of  b,c,e,  and  f; 
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mvMMm wmnMimm »!*»♦#  -r^j 


Ix-  -jk 


/■■(•• 

n- 


Aer ' 2+ ( Be-t-Af ) r ' +Bf 
(l~e2)r ' 2+2(b-ef )r'+c~f2 


dr  ' 


-Jk 


+ F, 


F1<rn+-rn-)  + V“  iTHi5? 

n-  n- 


- tan 


-1 


<l-e4)rn++(b-ef) 
l(l-e)  (c-fj-(b-ef)2 


O-e  )r  + ( b-ef ) 
n — 

^(l-e2(e-f2)-(b-ef)2 


(B-l 7) 


where 
F,  - D + 


Ae 

l-e‘ 


V - (Be+Af ) (l-e*)-2Ae(b-ef ) 

11 2 2 2 

2<l-eV 


F3  - [2Ae(b-ef)2~Ae(l-e2)(c-f2)-(Be+Af)(b-ef)(l-e2)+Bf(l-e2)2] 


* [(1-c2)2  ^1-e2) (c-f 2)-(b-ef )2]-1 

The  tabulated  Integrals  Dw  160.01,  160.11,  and  160.21  aid  in 
the  evaluation  of  1^1  I2  raay  be  evaluated  using  Dw . 380.001 
aa  r 


I2  - D(l+jkR 


■’/ 


n+ 


dr* 

R 


n- 


f 

The  abbreviations!^  and  OR  refer  to  Tables  of  Integrals  and 
other  Mathematical  Data.  Fourth  Ed.,  H.B.  Dwight,  Macmillan, 
N.Y.,  1961 ; and  Tables  of  Integrals,  Series  and  Products  I . S 
Gradshteyn  and  I.W.  Ryshik,  Academic  Press,  N.Y.,  1965,  res- 
pectively. 


7 5 


R ,+r  +b 

D(l  + JkRn)£n  ^"-±- 
n-  n- 


(B-18) 


The  subscripts  n,  n+,  and  n-  denote  quantities  which  are 
evaluated  at  r'-r^,  rn+,  and  r^_,  respectively. 

The  evaluation  of  I ^ is  facilitated  by  expanding  it 

in  partial  fractions  and  using  the  aubstitution 

J 2 

sinh  0 * (r’+b)/  c-b  to  obtain 


I3  - (1+jkR 


r n+ 

. if  SSMl 
J r2-c 


n- 


(1+jkR  ) 


) /*li+  r 

I — - 

J R"C 


_i_)  dr. 


R+C 


n- 


( 1+jkR 


n+ 


J ? 

(A  e-b4  sinh  0 + B-Ab) 


n - 


/ 2 / 2 
c-b  cosh  0 - e c-b  sinh  9 - f+be 


J 5 J 2 

c-b  cosh  0+  e c-b4  sinh  0 + f-be 


dO 


4 


Using  GR  2.451.2  and  GR  2,451.4,  one  finds  the  latter  integral 
to  be 
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I 


W!T9i  fWs* f Jftdi,vuwiiw<*m*> 


. R ,+r  ,+b 

1 3 - (l+JkRn)  ^e_  ln  iB±^L+-- 

1-e  n-  n- 


2 ( 1-e  ) 


(R  +-C  ,)(R„  +C  ) 
nr  nr  n-  n~ 

(R  -C  ) (R  ^+C  . ) 
n-  n-  n+  n+ 


+ tan 


, c-b2(R  -C  ) - (be-f ) R ,+b2-c 
- 1 n _n n+  

W+b> 


~tan 


. c-b  (R  -C  ) - (be- f ) R H>  -c 
■ 1 n-  n-  n- 

>5  <%.+») 


-tan 


-i  °~b  <R,++c,,+>+<t«-t>lt,,++b  -c 
F5(rn++k) 


+ tan 


, ^c-b2  (R  +C  )+(be-f)R  +b2-c 

- 1 n-  n- nj; 

P5(Cn-+b> 


F*  ' 


Hl-e  ) - Ae  1 be- i 
P 3 (1-e2 ) 


“ (c-b 2 ( 1-e 2 ) - (be- f ) 2 


Equations  (B-17 ) - (R-19)  complete  the  evaluation  of  (B-16) 
The  derivative  terms  appearing  In  (B-12)  require 
evaluation  of  integrals  of  the  form 


K dr*  - 


/■*(•*  T'?] 


( B — 2 0 ) 


where  the  unorltned  variable  r is  r^  or  r ^ , as  appropriate. 
The  two  edge  terns  also  have  non-integrable  singularities 
at  the  edge  which  cancel  between  the  various  terms.  To 
handle  this  situation  numerically,  the  singularity  must  be 
explicitly  identified  and  removed  for  numerical  integration. 
Thus  the  same  kind  of  approximate  analytical  integration  of 
(B-20)  as  used  to  evaluate  <B-16)  would  both  eliminate  one 
integration  and  explicitly  identify  the  singular  term.  The 
derivative  can  be  taken  inside  the  integral  if  care  is  taken 
to  Identify  the  oingular  terms.  Noting  that  8 C/3r  ■ 1, 

3 (R^-C^)/3r  ■ 0,  and  3A/3r  * 3B/3r  ■ 0,  we  have 


With  the  approximation  of  (B-15),  the  above  may  be  written 


-jkR 

e n(I4  + I3  + I6  + I?  + I8)  (B-21) 


I 

The  various  integrals  appearing  in  (3-21)  are  defined  and 
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Ipguw  jiiy  w «" 


j 


- k ‘ 


r 

/ 


n+ 


C ( er ' +f ) (Ar’+B) 


2 2 

R(R  -C  ) 


dr  ' 


n- 


i 

7 


n + 


.2  J CfAer  ' 2+(Af+Be) r '+Bf  1 , , 

“ ~ * 1 2 2 ar 
R(R  -C  ) 


n- 


Dividing  R2-C2  - (l-e2)r'2+  2 (b-e f ) r ' +c - f 2 into  the 
bracketed  term  in  the  numerator  of  the  integrand,  we  may 
write  the  integral  as 


n+ 


,2  f C I"  Ae  (Af+Be+W)r’+Bf+tf|  , , 

J » l-.2  R2-C2 


n- 


where 


w - _ 2Ae(b-ef) 

i 2 

1-e 


U - . Ae<;c_-L  l 

1 2 
1-e 


Expanding  the  second  term  in  brackets  in  the  Integrand  in 
partial  fxactions,  one  can  write  Ia  as 


I - I ' + I " 
» 8 8 


where 


V “ 


-k2Ae 

T 


1 -e 


i 

/ 


n+ 


e r 1 +f 


dr  1 


n- 


1 

1 

7? 

!• 

«<R  -R  ) + (f-beUn 

n*r  n - 

R +r  +b 
n+  n + 

1-e2 

R +r 

+b 

_ 

n- 

n- 

(B-26) 
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and 


X " . zX 
o 2 


i 


n+ 


[ (Af+Be+W)r  *+Bf+U] 


R-C  R+C 


dr_’ 
R ~ 


n- 


J o 

The  substitution  sinh  0 » (r'+  b)/  c-b  enables  one  to  write 


V a9 


T M m **R 
8 2 


'I 


n+ 


J 9 

[(Af+Be+W)  c-b  sinh  6 + Bf -HJ-b ( Af +Be+W)  ] 


n- 


J 1 J J 

c-b  cosh  0-e  c-b  sinhS-f+be 


J y ,/ y 

c-b  cosh  9+*  c-b  sinhS  + f-be 


dr  ’ 


Using  GR  2.451.2  and  GR  2.451.4,  we  obtain  finally 


v-  -k  < 


Af +Be+W 
2(l-e2) 


in 


( 1-e  ) 


R +c  +b 
n+  n+ 

R +r  +b 
n-  n- 


F 


7 U®11'1 


- tan 


,/<:-b2(Rn+-c„*>-<>>*-f)Rn++b2-c1 


^c-b2  (R  -C  ) - (be-f  ) R +b2-J 
n-  n-  «- 


(rn-+b)F5 


Can 


-I 


/^I(«„++c„+)  + l^-f)R„++b2-o 


<rn++b)F5 


+ tan 


-1 


c-b2(R  +C  ) 4-  (be-f)R  +b2-c 
n-  n - n- 


(r  +b ) F 
n-  j 


(B-27 ) 


where 


v - [Bf+U-b(Af+Be+W)1 ( 1-e  ) -e ( be- f ) (Af+Be+W) 

r i 2 

( 1- e ) F 5 

Equations  (B-22)  through  (B-27)  complete  the  evaluation  of 

the  integral,  (B-21).  Recall  that  the  Integral  (B-21)  needs 

to  be  evaluated  only  for  observation  points  at  the  bicone 

edge  (see  Eq.  (B-12)),  For  the  source  current  pulse 

associated  with  the  bicone  edge,  there  results  a non- 

integrable  singularity  (with  respect  to  <ji  ’ integration) 

which  comes  from  the  term  l/R  . in  (B-24).  Each  of  the 

n+ 

derivative  terms  in  (B-12)  contains  such  a non-integrable 
singularity,  however,  and  they  are  of  opposite  signs  so  as 
to  cancel  each  other.  For  numerical  integration,  of  course, 
the  canceling  singularities  must  be  analytically  subtracted. 

The  integrals  1^  through  tg  contain  integrable  singu- 
larities such  as  the  usual  one  where  source  and  field  points 
coincide  (i.e,,  R-0) . In  addition,  however,  there  are  also 
Integrable  singularities  introduced  by  the  transformation 
(B-6).  These  arise  from  current  sources  which  lie  along  and 
are  directed  transverse  to  the  line  which  passes  through 
the  observation  point  and  which  is  in  the  direction  of  the 
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electric  field  component  of  interest.  For  example,  in 

computing  the  tangential  field  component  F.  on  the  cone, 

c 

singularities  arise  from  currents  on  the  topcap  and  its 
image.  While  the  many  singularities  complicate  the  numerical 
procedure,  they  can,  in  principle,  be  handled.  However, 
the  unwieldiness  of  the  functions  appearing  in  the  integrals 
I.  through  IQ  makes  the  numerical  procedure  rather  ineffi- 
cient  and  subject  to  error.  For  example,  the  computer  code 
derived  from  the  formulation  presented  here  seemed  to  yield 
reasonable  results  for  moderate  cone  angles,  but  often 
yielded  erroneous  results  for  the  very  small  cone  angles 
used  to  check  the  program.  The  complexity  of  the  formulation 
became  a considerable  hinderance  in  determining  the  source 
of  these  difficulties.  Consequently,  the  final  calculations 
were  done  using  the  formulation  of  Section  III  which  was 
developed  as  an  extension  of  methods  currently  being  used 
to  treat  flat  plate  surfaces. 
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